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One of the most successful results of perturbative QCD is surely the detailed de- 
scription of logarithmic scaling violations in Deep Inelastic Scattering (DIS). Recent im- 
provements in experimental accuracy are now opening up the possibility of measuring 
subleading perturbative corrections, as well as corrections to perturbation theory as ex- 
pressed by higher twist effects. However, it was pointed out in that in the region of 
~ 1 (GeV/c)^ the scale dependence of the non-singlet nucleon structure function may 
be qualitatively rather different from that predicted by purely perturbative QCD, due to 
non-perturbative effects related to the anomalous breaking of axial U(l) symmetry. 

In the effects on non-singlet structure functions evolution were studied, producing 
testable predictions of this model. This predictions are in good qualitative agreement with 
the experimental data available today and could lead to a clear confirmation of the model, 
when more precise data will be available. 

Here we will study the effects on the singlet structure function evolution, which are 
expected in this scheme 0. 

Let us first summarize the main features of the formalism [0. 

In the model of 0] is introduced, beyond the usual Altarelli-Parisi QCD evolu- 
tion, the radiation of bound states (pseudoscalar mesons) by quarks. This leads to modify 
the evolution equations, adding new splitting functions, functions that describe this phe- 
nomenon. 

In Ref. the mesons (H(/c)) are coupled to quarks by the effective action 

where U{k) = exp (z75H(A;)//n), /n is the meson decay constant and the coupling is 
defined by 

X{k,p) ^ S-\p + ^k){Omp - \k)i^{p + \k)\Il{k))S-\p - ifc), (2) 
where S~^{j>) = + is the quark propagator (S(p^) is the self-energy). 



^ Incidentally, a previous tentative qualitative estimate of these non-perturbative effects on 
singlet evolution from the results obtained in the non-singlet case had led to the suggestion |^ 
that the effect could be opposite in sign and larger than the one we calculate here, giving a possible 
explanation to the albeit marginal difference between the Os extracted by DIS and LEP data Q. 
This quantitative analysis does not confirm the earlier qualitative estimate. 
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The effective couphng xi^^p) is then expanded, in the most general case, in terms of 
four vertex functions (a pseudoscalar, two axials and a tensor.) 

This permits us to evaluate the cross section a'^.^ {x; t) for the absorption of a virtual 
photon 7* and the emission of a pseudoscalar meson X and, then, the splitting function 



For the sake of simplicity, considering that the mixing angle is small (see, for example, 
0), we assume in the following that r] and r]' are the pure octet and singlet respectively. 
One obtains for the non-diagonal {V^q^ , flavour changing) and for the diagonal (^^) 
splitting functions: 

^.T = -Pi"" = ^a^*-^ = ^a^*^~ 
""^ '^'^ dt dt 



The strange mesons [K) contributions are kinematically suppressed and therefore ne- 
glected. In the following we will thus neglect anomalous contributions to the strange 
sea evolution. 

In the distribution 

q^{x) = u{x) + u{x) — d{x) — d{x) (5) 

was studied. 

For this the (mesonic) non-singlet evolution is given by 



d 

di' 



u-d)+ = {Vg-V^q^)®{u-d)+. (6) 



Because of the strong mass difference between rj' (m^/ = 958 MeV/c^) and pions 
(rriT^o = 135 MeV/c'^), which is caused by U{Nf) breaking due to axial anomaly 0, there 
is a relevant difference between "P^ and P^^, which is responsible of quark-antiquark sea 
SU{2)p symmetry breaking. 

The fact that the U{Nf) breaking concerns pseudoscalar mesons only implies that 
these bound states only contribute to non-singlet evolution: for other mesons' multiplets, 
V^g and are almost equal and moreover they are kinematically suppressed (of about 

one order of magnitude) due to the larger masses of these mesons. 
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The explicit calculation of cross sections a'^.^{x;t) |l|] shows that they depend on 
scalar (99, which dominates at intermediate Q^) and axial {(p, which determines the 
tail) vertex functions only. 

A dipole form is then assumed for the vertex functions {p is the quark momentum, 
U = 93 MeV) 

^(P^) = 'f^. (8) 

in terms of the three parameters A, A and (7^. In A = A is hypothesized (the dependence 
of the results on the precise value of A is small anyway), while (7^ is between and 1. 

The A value may be related, using chiral Ward identities, to the constituent quark 
mass. This requires it to vary in the interval 0.4 GeV/c^ ^ A ^ 0.8 GeV/c^. The A and 
Qt^ values have been fixed in Refs. |]l|J^ to A = 0.4GeV/c^ and g.,^ = 0.5 by the request of 
reproducing the NMC Gottfried sum datum 0. 

In |jl|] it has also been calculated that the anomalous dimension due to pion emission 
is significantly different from zero in the region 0.05 (GeV/c)^ ^ ^ 5-10 (GeV/c)^. 
The evolution must therefore be extended up to very low Q^. 

The strategy of was to insert the effect of the meson emission inside the model of 
Ref. [0, which permits an extension to the evolution up to very low in order to 
study the non-singlet evolution starting from a reasonable ansatz at ~ 0(GeV/c)^. 

In the following we will adopt a different strategy in order to study the anomalous 
singlet evolution, namely we will begin the evolution at an input scale Qq sufficiently 
high to allow other non-perturbative effects on evolution to be neglected a part from 
the meson emission ones (in the following we will assume as a reasonable choice ^1- 
1.5 (GeV/c)^). This allows us also to use as input an ordinary parametrization of the 
experimental structure functions. 

As we are evaluating a first approximation of these effects, we will neglect (as in [|l],H) 
the effects due to mesons other than the pseudoscalar ones, which are expected to be 
kinematically suppressed due to larger masses. 



For the singlet evolution one has to consider the full set of coupled equations (omitting 
gluon contributions for simplicity): 

j j a 



dt 

d 
di 



where n"(x;t) are the distribution functions of mesons of the pseudoscalar multiplet. 
The distributions must fulfill the momentum sum rule 



1 

dx ■ X 







(10) 



where g is the gluon distribution. 

The first two equations in (^ express the evolution of the quark and antiquark dis- 
tribution due to mesons emission, whereas the last equation gives the evolution of the 
pseudoscalar meson distribution. The former are determined by the splitting functions 
Vqq, which expresses the probability of a quark to be emitted by another quark (with 
emission of an unobserved meson), and "Pgn, which expresses the probability of a pseu- 
doscalar to fragment into a quark and an antiquark (one of which is observed) ; at first order 
Vqq is vanishing. The latter is determined by the splitting functions Vnq and Pnn, which 
give the probability of a pseudoscalar to be emitted by a quark or another pseudoscalar, 
respectively. 

It must be noticed that not all these splitting functions are independent, being related 
by isospin and charge-conjugation invariance (see |1[]). 

Considering that we will be interested in the medium-to-large x region (far from the 
low X region where one expects a rapid growth of the singlet structure function due to the 
radiative generation of the sea from the glue), we will neglect, as a first approximation, the 
effect of the splitting of the radiated mesons in a qq pair, which is expected to contribute 
substantially to the low x region only and that is in any case also suppressed by the fact 
that the mesons distributions are not expected to be very large [||. 



We will impose the momentum conservation, adding, as usual [10|, a delta contribution 
to the splitting functions: 



p^(^)=p^(^)-p^(5(z-l) (11) 

and 

P^^{z) =V^^iz) -VaS{z-1), (12) 

where ^ 

Va= [ dz-z-{V''{z)+V''''{z)). (13) 

The fact that we have neglected the mesonic distributions is thus phenomenologically 
compensated, in order to keep the validity of the momentum sum rule (p!0|), by the explicit 
choice of the value of Pa- 

The F2 evolution is then given by: 

= [ J [PApix/y) + P^{x/y) + P^^(x/,)] X 

[4/9 ■ (uiy) + u{y)) + 1/9 • {d{y) + d{y))\ + ^l^) 
PAp{x/y)xl/%-[s{y) + -s{y)], 

where the mesonic effects concern the u and d quarks only (as we are neglecting the 
kinematically suppressed contributions of strange mesons) and Pap ai's the usual Altarelli- 
Parisi splitting functions ^ . The heavy quarks, considering the low under investigation, 
are completely negligible. 

The numerical calculations are made starting with an input for the different parton 
distributions at Qq = 1.6(GeV/c)'^; in the following we will show the results obtained 



using as input the CTQE3M distributions ||TlI]; however, we have explicitly checked that 



our results are largely independent from this choice. 

We carry out the evolution (|l^) through many small steps in (we have explicitly 
verified the stability of our results on the choice of the steps), adding at every little step in 
(5^ the mesonic contribution to the one coming from the Altarelli-Parisi part (at leading 
order). This last part of the evolution is effected in the momentum space (in order to have 
a better numerical accuracy); therefore, at each step, the singlet, non-singlets and glue 
distributions are parametrized in the form pi -x^^ . (^\ — xY^ ■ {l+p^x'^^^ -{-p^x+p^x^ +p-i:x^), 
which permits an analytical Mellin transform to the momentum space. The inverse Mellin 
transform to x space is made numerically. The value A^q""^ = 200 MeV/c is used. 
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The mesonic effects tend to increase the growth of the parton distributions at low x 
(under x ~ 0.15), while at larger x the effect of the delta term dominates producing a 
more negative '^^'^^^^ ■ Altogether this therefore simulates a larger value of ^^(Q^), except 
for a very narrow region around the point where the mesonic contribution changes of sign. 

In fig. |l| '^^"^^'^^ is shown, as obtained with (solid curve) and without (dashed curve) 
mesonic effects at = 1.6(GeV/c)2. 

For X ^ 0.6-0.7 there is a difference of 0.03 (~ 10%) between the two predictions (at 
higher x non-leading twist effects, as target mass corrections, may be quite important). At 
this scale the mesonic effects simulate an increase of A2o"^ from 200 MeV/c to 210-220 
MeV/c. At larger scales the effect reduces rapidly. It is already small at ~ 5 (GeV/c)^ 
(where it is reduced by about one half of its former value), albeit always leading to an 
overestimate of Aqcd of ~ 5%, and it is negligible for ~ 10 (GeV/c)^. 

It is thus not yet really possible to observe this effect in today's DIS data; anyway the 
order of magnitude of the effect at ~ 1 (GeV/c)'^ is such that it could constitute a non 
completely negligible source of error in the determination of by DIS data. A precise 
determination of this error obviously depends on the weights assumed by different and 
X regions in a global fit; however, one can roughly expect a possible overestimate of Aqcd 
of a ^ 5-10% for a global fit in a region such as 1 (GeV/c)^ ^ ^ 5 (GeV/c)^ and 
X ^ 0.01. 

When higher precision experimental data will become available it would surely be 
quite interesting to search for a confirmation of this effect by comparing '^'"^'^^ evaluated 
with Aqcd obtained in other processes and the one extracted by DIS data in the ^ 
5(GeV/c)^ region. However, a direct observation of this effect could turn out to be rather 
difficult due to the fact that scale fixing and other theoretical uncertainties could lead to 
an arbitrariness in determining a precise Aqcd value, which is of the same order as this 
effect. 

The study of the singlet anomalous evolution enables us to evaluate the anomalous 
evolution effects on "'^^^^(^^ also (where F2 is the neutron structure function and F2 the 
proton one.) 

It must be noticed that the numerical estimate '^^^^J-^^ is much more sensitive to the 
choice of the parton distributions than the evaluation of the derivative of F2. The estimate 
of the effect of the mesonic contributions to this quantity is thus, for a precise numerical 
calculation, partially depending on this choice as well. Anyhow for a first evaluation of the 
effect this dependence is not really severe, as we have checked with other inputs beyond the 
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CTEQ3M parton distributions (used as input at = 1.6(GeV/c)^ as in the derivative 
case). 

In this case the mesonic effect is a positive contribution to '^'^^^l^^ ^ in the whole range 
of X. It enhances this quantity with respect to the one evaluated with only the Altarelli- 
Parisi splitting functions of ~ 0.015 at = 1.6(GeV/c)^ in the region 0.1 < x < 0.6. The 
effect is still larger at higher x, where a careful analysis of target mass effects and of other 
higher twist effects would anyway be necessary, while it disappears at lower x (at x = 0.01 
it is reduced to a difference of 0.003). 

The effect decreases with and it is reduced by 50% at = 10(GeV/c)^. 

Nowadays experimental data [|12| are not yet able to discern this effect because of 
large errors. Anyway they are in qualitative agreement with this prediction for x ~ 0.4- 
0.6, where they show values larger than zero, above QCD predictions, while they seem to 
point out more negative values than the QCD prediction for x ~ 0.1-0.3. In |jl^ this effect 
has been interpreted in terms of higher twist effects. 

More precise experimental data and a more careful treatment of higher twist effects 
(for example by fitting them carefully in high precision experimental data or by evaluating 
them in some theoretical scheme) could lead in a next future to an at least qualitative test 
of this prediction. 

In summary we have presented an estimate of the effects of the anomalous evolution 
due to the mesonic degrees of freedom on singlet structure functions in the scheme of |jl[]. 
The effects have been found to be not very large; nevertheless they are not negligible for 
1 (GeV/c)^ ^ Q2 ^ 5 (GeV/c)^ and could constitute a non-negligible source of error in a 
high precision determination of from DIS including data in this region. 
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Figure Captions 

Fig. 1. The ratio ^IoIq^ obtained with (sohd curve) and without (dashed curve) mesonic 



effects at Q'^ = 1.6(GeV/c)^. The value Al{^^ = 200 MeV/c is used 
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